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The persistent current in a lattice model of a one-dimensional interacting electron system is sys- 
tematically studied using a complex version of the density matrix renormalization group algorithm 
and the functional renormalization group method. We mainly focus on the situation where a sin- 
gle impurity is included in the ring penetrated by a magnetic flux. Due to the interplay of the 
electron-electron interaction and the impurity the persistent current in a system of N lattice sites 
vanishes faster then 1 /N. Only for very large systems and large impurities our results are consistent 
with the bosonization prediction obtained for an effective field theory. The results from the density 
matrix renormalization group and the functional renormalization group agree well for interactions 
as large as the band width, even though as an approximation in the latter method the flow of the 
two-particle vertex is neglected. This confirms that the functional renormalization group method is 
a very powerful tool to investigate correlated electron systems. The method will become very useful 
for the theoretical description of the electronic properties of small conducting ring molecules. 
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I. INTRODUCTION 

The experimental observation of persistent currents in 
mesoscopic metallic and semiconducting rings pierced by 
a magnetic flux 1-6 has led to many theoretical inves- 
tigations focusing on the interplay of electron-electron 
interaction and disorder in such systems. 7 This inter- 
play is considered to be one of the possible reasons 
for the large current observed in the experiments. 7 De- 
spite these studies a quantitative theoretical understand- 
ing of the observed amplitude of the currents for the 
three-dimensional rings is still missing. To gain theo- 
retical insight the simplified situation of one-dimensional 
rings with interaction and disorder was studied using 
exact diagonalization for systems of very few lattice 
sites (up to 16) 8-10 and the self-consistent Hartree- 
Fock approximation. 10,11 Here we consider the further 
simplified problem of the persistent current in a one- 
dimensional ring of interacting electrons in the presence 
of a single impurity and penetrated by a magnetic flux. 
Within an effective continuum field theory it has theoret- 
ically been investigated using bosonization 12 and confor- 
mal field theory 13,14 At first glance this problem seems 
to be of purely academic interest, but the fast progress in 
the design and manipulation of conducting ring molecules 
gives a perspective that such systems might be accessible 
to experiments in the very near future. 

For a continuum model of non-interacting one- 
dimensional electrons the leading behavior of the persis- 
tent current I{<ft) in the system size L can be calculated 
in the presence of an arbitrary potential scatterer. 12 It 
is a periodic function of the flux, vanishes as 1/L, and 
its shape and size are determined by the absolute value 



of the transmission amplitude |T(fc^)| of the potential 
at the Fermi wave vector fci?. For a vanishing impurity, 
i.e. |T(fcj?)| — > 1, I {<ft) has a saw tooth like shape, which 
gets rounded off if \T{Uf)\ is decreased. In the limit of 
small |T(fcp)|, I {(ft) is proportional to |T(&f)| sin.0. For 
the tight-binding lattice model supplemented by a single 
weak hopping matrix element the persistent current at 
half-filling can also be calculated analytically 13 and the 
same characteristics can be found (see Sect. II). 

Compared to the Fermi liquid behavior of higher di- 
mensional systems a large class of models of homogeneous 
one-dimensional interacting electrons has a significantly 
different low-energy physics. These models belong to the 
Luttinger liquid universality class. 15 The low-energy ex- 
citations of Luttinger liquids are not given by fermionic 
quasi-particles, but are of collective, bosonic nature. This 
leads e.g. to a typical power-law decay of correlation 
functions. The low-energy physics of Luttinger liquids is 
characterized by a set of interaction and filling factor de- 
pendent parameters. 15 For the case of spinless fermions 
on which we focus any pair of the four parameters vj 
(velocity of current excitations), vn (velocity relevant if 
particles are added) , v c (velocity of charge excitations at 
constant number of particles), and the Luttinger liquid 
parameter K can be used. Within bosonization and for 
the impurity free case the persistent current in a Lut- 
tinger liquid is periodic in <\> and of saw tooth like shape 
with slope vj/{ttL). 15-17 In Sect. Ill we will compare our 
numerical results to this prediction. 

The low-energy physics of Luttinger liquids is strongly 
affected by the presence of a single impurity. 18-23 The 
problem is usually mapped onto an effective continuum 
field theory using bosonization, where terms which are 
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expected to be irrelevant in the low-energy limit are 
neglected. 18-22 Within this field theory the leading L de- 
pendence of the persistent current was obtained by an 
additional self-consistent approximation using the anal- 
ogy to the problem of quantum coherence in a dissipative 
environment. 12 This approach gives a current which for 
large L vanishes as L~ aB ~ x , with as = l/K — 1, and 
independent of the bare impurity strength is of purely si- 
nusoidal shape. 12 For repulsive interaction one has K < 1 
and thus as > 0. as is also the exponent of the power- 
law suppression (as a function of energy) of the local spec- 
tral weight close to an open boundary and the chemical 
potential. 21,24 This explains the index B which stands 
for boundary. Many of the bosonization results for ob- 
servables which are dominated by the interplay of a sin- 
gle impurity and the electron-electron interaction can be 
understood in terms of a single particle picture, which 
for Luttingcr liquids has of course to be used with cau- 
tion: For large L the effective transmission amplitude 
near the Fermi points is suppressed with respect to the 
non-interacting transmission by a factor of L~ aB . Com- 
bining this with the behavior of the persistent current in 
the non-interacting case the above result, obtained with- 
out the use of the single particle language, can be derived. 

For several reasons it is desirable to directly show 
the above behavior of the persistent current in micro- 
scopic lattice models avoiding bosonization. Mapping 
such models on the field theory involves approximations. 
Their validity can be questioned and they lead to a loss 
of information; the scales of the microscopic model and 
corrections to the expected power-law scaling hidden in 
the irrelevant terms are lost. Furthermore even within 
the field theory an additional approximation is necessary 
to determine the leading (in the system size) behavior 
of the persistent current. A knowledge of the detailed 
shape of I(<p) beyond the leading behavior and for mi- 
croscopic models is of special importance if one wants to 
compare theoretical results to experiments. Several at- 
tempts have been made in this direction using the model 
of spinless fermions with nearest neighbor hopping and 
interaction. Here we also focus on this model. The persis- 
tent current can be calculated by taking the derivative of 
the groundstate energy E (4>) with respect to the flux 4> 
penetrating the ring. Instead of calculating the full func- 
tional form of I (<j>) the so called phase sensitivity AE , 
which is the difference of the groundstate energy at flux 
and 7T has numerically been determined using the den- 
sity matrix renormalization group (DMRG) method. 25 ' 26 
The phase sensitivity can be considered a crude measure 
for the persistent current, but of course does not con- 
tain information on the detailed shape of the current as 
a function of </>. Very recently AE has again been stud- 
ied numerically using DMRG and the quantum Monte 
Carlo method. 27 AE instead of /(</>) (which implies cal- 
culating E (<j)) for several </> and numerically taking the 
derivative) is calculated because this way the hamilto- 
nian matrix remains real and the numerical effort is re- 
duced considerably. 28 Here we use complex DMRG to 



calculate I{<j)) with very high precision and for systems 
of up to N = 128 lattice sites. 29 Additionally the func- 
tional renormalization group (RG) method introduced re- 
cently into the theory of strongly correlated electrons is 
used. 30 ' 31 It has been applied to two-dimensional corre- 
lated electron systems 32 and one-dimensional Luttingcr 
liquids. 33-35 We have used this method before to study a 
local observable 33 ' 34 - the local spectral weight close an 
impurity - which is also dominated by the interplay of 
electron-electron interaction and impurity. In contrast 
the persistent current is a property of the entire sys- 
tem. Within the RG approach the flow equations are 
closed by neglecting the flow of the two-particle vertex. 
Nonetheless DMRG and RG agree quantitatively for in- 
teractions of the order of the band width. This con- 
firms that the functional RG is a very powerful tool to 
investigate strongly correlated electrons. For the single 
impurity case we indeed find that the persistent current 
vanishes faster then TV -1 . To analyze NI(<t>) in more 
detail we expand the current in a Fourier series, demon- 
strate that for large impurities and very large system 
sizes the behavior of the first Fourier coefficient is con- 
sistent with a power-law decay with exponent —as, and 
that the higher order coefficients decay even faster. Thus 
in the ./V — > oo limit I(4>) is proportional to sin^ with 
an amplitude which vanishes as N~°' B ~ 1 . However, this 
universal bosonization prediction only holds for very long 
chains respectively very large impurities. For smaller sys- 
tems and impurities the asymptotic limit is not reached 
and the current displays a more complex behavior. It 
can quantitavely be described using the functional RG 
method. This makes the method an ideal tool to inves- 
tigate the electronic properties of small, one-dimensional 
molecular rings which might be experimentally accessible 
very soon. 

II. THE MODEL 

The Hamiltonian for the impurity free ring penetrated 
by a magnetic flux <f> (measured in units of the flux quan- 
tum <po — hc/e) with nearest neighbor interaction U is 
given by 

H = -±(c]c j+1 e^ + c] +lCj e-^) 

JV 

+ u ^2 n i n i+i: (!) 

in standard second-quantized notation. The hopping ma- 
trix element and the lattice constant are set to one and 
periodic boundary conditions are used. For U = the 
groundstate energy of the model can easily be calculated. 
At temperature T = and for fixed N the persistent 
current follows from Eo(cj)) by taking the derivative with 
respect to <j) 
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as can be shown using the Hellman-Feynman theorem. 
To leading order in 1/N this gives 

,, ... _ Vf J for A^podd and — tt < <fi < tt 
W ~~ ~t7n X \(j)-TT for N F even andO < <f) < 2tt , 

where vp denotes the Fermi velocity and Np is the num- 
ber of particles. 36 Both functions have to be continued 
periodically. Eq. (2) also holds for non-vanishing inter- 
action and if impurity terms are added. The above even- 
odd effect can also be observed if interaction and impu- 
rities are included. We from now on only considere even 
N F . 

At = the Luttinger liquid parameter K and the 
velocity of current excitations vj of the model Eq. (1) 
can be determined using the Bethe ansatz. 37 For half- 
filling the resulting integral equations have been solved 
analytically with the results 



K = 



and (vj — v c K) 
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arccos (-U/2) 



Vj = TT 



y/i - mm 2 

tt — arccos (-U/2) 



• arccos (— U/2) 
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For this filling the model is a Luttinger liquid for — 2 < 
U < 2. At \U\ = 2 the model shows phase transitions 
to a charge density wave groundstate (U = 2) and a 
phase separated state (U = —2). To leading order in 
1/A the bosonization prediction for the persistent cur- 
rent in a homogeneous Luttinger liquid is a saw tooth 
like curve with slope —vj/(ttN), i.e. compared to the 
non-interacting case vf is replaced by the velocity de- 
termining the current excitations. To compare our data 
with this result also away from half-filling we numerically 
solved the Bethe ansatz integral equations following Refs. 
37 and 38 and determined vj (see Sect. III). In our study 
we will always stay in the Luttinger liquid phase. 
To if we add a hopping impurity 



H h = (l- p) {cW iHN + c\c N e-^ N ) . 



(5) 



with p between (no hopping between sites N and 1) 
and 1 (no impurity). In the non- interacting limit, at half- 
filling, and at wavevector Uf the transmission coefficient 
of the hopping impurity (for N — ► oo) is given by 39 



|T(MI 2 



,2^2 
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which provides us with a measure for the strength of the 
impurity. For U = 0, Np = N/2, and to leading order 
ml/N the persistent current of the hamiltonian H + Hh 
has been calculated using conformal field theory. 13 To 



obtain I((f>) it is not necessary to use this field theoretical 
approach since Eq. (7) of Rcf. 13 which determines the 
allowed wave vectors can be solved directly. Written in 
terms of \T(kp)\ the current is given by 

m = ^ arCC ° S (|T(MI C ° S [<t> " " ]) \T(k F )\ sin0 , (7) 
V 7 ttN ^i- |T(M|2 C os 2 

which is exactly the expression obtained in Ref. 12 for 
the continuum model. Later we will be interested in the 
small T(fcp) limit. Therefore we expand Eq. (7) up to 
third order in |T(fc F )| 

NIW - V -j\T(k F )\ (1 + |T(fc F )| 2 /8) sin(0) 
+g|T(M| 2 sin(20) 

+^|T(M| 3 sin(30) + O(|T(fc F )| 4 ) , (8) 

which at the same time gives an expansion in a Fourier 
series. The expansion explicitly shows that in the limit 
|T(fcp)| — > of a strong impurity the current becomes 
more and more of sinusoidal shape. 

A simple approximation which allows to study the ef- 
fect of the interaction and impurity on the persistent cur- 
rent simultaneously is the Hartree-Fock approximation. 
For the bulk properties of homogeneous one-dimensional 
correlated electron systems this approximation does not 
capture any of the Luttinger liquid features and is thus 
of very limited usefulness. Furthermore, when apply- 
ing the self-consistent Hartree-Fock approximation to a 
model with a single impurity, the self-consistent itera- 
tive solution of the Hartree-Fock equations will drive the 
system into a charge density wave groundstate with a fi- 
nite single-particle gap which is qualitatively incorrect 
since a single impurity cannot change bulk properties 
of the system. Nevertheless self-consistent Hartree-Fock 
has been used to determine the persistent current in one- 
dimensional, interacting, and disordered rings 10 ' 11 and 
also for the single impurity case. 11 In contrast if one is 
interested in the local properties close to a boundary or 
impurity non-sclf-consistent Hartree-Fock provides use- 
ful informations. 24 > 33 > 34 In Sect. V we also present re- 
sults for the persistent current calculated within the non- 
sclf-consistent Hartree-Fock approximation. They were 
obtained by numerically determining the groundstate of 
H + Hh for U = and fixed <j> and N. From this the ex- 
pectation values (rij) and ^c] +1 Cj^ , which determine 

the mean- field hamiltonian Hmf , can be calculated. Hmf 
can then be diagonalizcd numerically. There are two pos- 
sibilities to determine the groundstate energy within this 
approximation. One can either take the Slater deter- 
minant expectation value of if + Hh using the ground- 
state of Hmf or determine the energy via the one-particle 
propagator 40 by the formula 
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MF 
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3 = 1 
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-E(^"<^«>„ + -) 
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(l-p)( e ^(4 Cl ) MF+ c.c.) 



(9) 



where the expectation values are taken using the ground- 
state of Hmf- Cj-^(r) denotes the annihilation (creation) 
operator in the imaginary time Heiscnberg representa- 
tion. Transforming the first term into the Matsubara 



frequency representation it can be written as £ i 
with the eigenenergies £^ F of -Hmf- Only if the self- 
consistent Hartree-Fock approximation is used, which for 
the reason given above we do not consider in the present 
context, both possibilities give the same result. Formally 
both approximations to the energy are correct to leading 
order in U , but it turned out that the latter method gives 
better results compared to the high-precision DMRG 
data. For the results presented in Sec. V we thus used 
Eq. (9). There it will be shown that the currents cal- 
culated using the Hartree-Fock approximation are qual- 
itatively wrong, since in this method correlation effects 
are neglected which in the present context arc of great 
importance. 



-MF 



III. COMPLEX DMRG 

Using the DMRG algorithm the groundstate energy of 
an interacting one-dimensional many fcrmion system can 
be calculated to high precision. 41 To determine the per- 
sistent current using Eq. (2) for the hamiltonian H + Hh 
[Eqs. (1) and (5)] the DMRG procedure has to be gen- 
eralized to complex hamiltonian matrices. Calculation 
time scales up by about a factor of 4, memory usage by 
a factor of 2. This limits the performance of the method, 
which is however numerically very stable as hermiticity 
is conserved. We have kept up to 400 states, ensuring 
that energies and derived currents are essentially exact. 

As a test of our program we first studied the impurity 
free case given by the hamiltonian Eq. (1). We calculated 
the groundstate energy as a function of <f> for < cf> < it. 
Results for quarter-filling n = 1/4, TV = 64 and U = 
0.5, 1, 1.5 are shown in Fig. 1. Bosonization predicts that 
to leading order in 1/N the current is given by 



/(0=-^(*-T) 

and thus the groundstate energy by 

E (4>) = const. + -^{4> - tt) 



(10) 



(11) 



As shown in Fig. 1 the DMRG data nicely lie on curves 
of this form with the constant and the current velocity 



as fit parameters. In Fig. 2 the v® M (n, U) extracted 
from such fits are compared to the exact vj(n,U) ob- 
tained from the Bethe ansatz. Data for n = 1/4 (with 
N = 64) and n = 1/2 (with N = 60) and different U are 
shown. For both fillings the Bethe ansatz and DMRG 
results are indistinguishable on the scale of the plot. The 
relative error of the DMRG data is of the order of 10~ 4 . 
This confirms the bosonization prediction, shows that the 
complex DMRG can be applied successfully, and that 
the current velocity can to very high precision be ex- 
tracted from finite size data as small as 60 lattice sites 
without using finite size scaling. In Rcf. 25 data for the 
phase sensitivity of the translational invariant model at 
half-filling were obtained using DMRG. Our results for 
iij MRG (n = 1/2, U) are consistent with the ones pre- 
sented there. As seen in Fig. 2 the dependence of vj on 
U gets weaker if the filling is getting smaller. This hap- 
pens because for smaller fillings the lattice model is closer 
to the electron gas model with quadratic dispersion. The 
latter model is Galilean-invariant which leads to vj = vp 
independent of the interaction. 
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FIG. 1. Groundstate energy as a function of the flux cf> for 
quarter-filling, N = 64, and different U. The symbols are 
DMRG data and the lines are quadratic fits (see the text) . 
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FIG. 2. Current velocity vj(n,U) 
Bethe ansatz and from the DMRG. 
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IV. RG METHOD 

Besides DMRG we will also use the functional RG 
method in the version using one-particle irreducible ver- 
tex functions. 30 ' 31 In collaboration with W. Metzner 
and K. Schonhammer we have successfully applied this 
method in the past to determine the local spectral weight 
of a Luttinger liquid close to an open boundary and an 
impurity. 33 ' 34 In the method one introduces a cut-off pa- 
rameter A in the non-interacting propagator G° cutting 
out degrees of freedom on energy scales less than A and 
derives an exact hierarchy of coupled differential flow 
equations for the one-particle irreducible vertex functions 
by differentiating with respect to A, where A flows from 
oo to 0. 

For <j> = the flow equation for the selfenergy 42 of the 
hamiltonian H + Hh has been given in Ref. 34. As in this 
reference we here also neglect the flow of the two-particle 
vertex which closes the set of differential equations; this 
leads to a energy independent selfenergy Within the 
approximation the results obtained are at least correct 
to leading order in U, but our work presented in Refs. 
33 and 34 shows that the Luttinger liquid scaling of the 
impurity (i.e. the transmission) is included, which makes 
the RG a promising method also in the present context. 

If a frequency cutoff 



g°< a m = e(M - a)g°m 



(12) 



is used, the set of flow equations for the selfenergy 34 in 
the Wannier basis is given by 



dA 



« V A 
dA^ 1 



2ir 



s=±lw=±A 



u 

2^ 



w=±A 



with 



G A (iuj) = | [G°(iw)] 1 -S A }" 



(13) 
(14) 

(15) 



In Eqs. (13) and (14) the flow from A = oo down to 
a scale A much larger then the band width has already 
been included. The flow is continued from Ao downwards 
with the initial conditions 



U 



^3,3 
^3,3±1 - U 



1 < j < N 
- ■ 1<3<* 

s tv = (i - p) ^ lHN = 



(16) 



In contrast to the situation studied in Refs. 33 and 34 
the matrix elements £< A j±i are complex numbers which 
increases the size of the set of differential equations. Fur- 
thermore we here have to use periodic boundary condi- 
tions so that G A (iuj) is not tridiagonal. 



The set of Eqs. (13) and (14) is complemented by a dif- 
ferential equation for the "zero-particle" vertex 70 (with 
the Boltzmann constant ks set to one) 



d 

dA 



lim T^ 7o A 



^Tr 



53 M 1 -G°(iuj)Z A (iLu)} 



Lw = ±A 



. (17) 



where Tr denotes the trace over the indices of the Wan- 
nier basis states. The initial condition is Jq = °° = Jq" = 
0. Eq. (17) also holds if the flow of the two-particle ver- 
tex and all higher vertices are taken into account. The 
flow of 7 A only couples to the selfenergy. It does not feed 
back to any higher order vertex function, 
related to the grand canonical potential ft by 



7o = 7o A= ° is 



n = t 7o + n 



(18) 



where fio is the grand canonical potential at U = 0. From 
f2 the ground state energy can be obtained in the T — > 
limit 



Eo(4>) =fi(T = 0) 



lhnT7o + £ o °(0)-Mo (n) q 



(19) 



A" is the particle number operator and \x the chemical 
potential. The quantities with an additional index are 
taken at U = 0. Eqs. (18) and (19) make explicit that the 
functional RG is a grand canonical method. In general 
to obtain the persistent current within a grand canonical 
ensemble corrections to Eq. (2) have to be taken into ac- 
count since the chemical potential can dependent on the 
flux. 43,44 For half-filling and a single hopping impurity \x 
is always fixed at U (and ^ at 0) and thus we do not 
have to worry about this issue if we use the RG only for 
n = 1/2. 

The set of differential equations (13), (14), and (17) 
can be solved numerically for a fixed and <j>. Addition- 
ally Eq(4>) can be calculated by numerically diagonalizing 
the one-particle problem at U = 0. Thus Eo(<p) can be 
determined. There are two factors which limit the system 
sizes that can be treated numerically. On the right hand 
side of the flow equations a, N x N matrix has either to 
be inverted or to be diagonalizcd. This limits JV to a few 
thousand lattice sites. The more restrictive limit comes 
from the fact that within the ground state energy which 
is proportional to A^ one is interested in the <f> depen- 
dence which vanishes faster than 1/N (see above and the 
next section). Thus the differential equations have to be 
solved with extremely high precision. For this reason we 
here only consider systems of up to A^ = 256 lattice sites. 
It is worth noting that for a given A^ the RG method 
with the approximation as we use it is still much faster 
than the "numerically exact" DMRG. 

An alternative way to determine Eq (<j>) would be to use 
Eq. (9) with the mean-field expectation values replaced 
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by the ones determined at A = 0. Formally the use of the 
flow equation for 70 and this method are correct to lead- 
ing order in U but it turned out that calculating E (<j>) 
from 7^ =0 gives much better results compared to the ac- 
curate DMRG data. For the results presented in the next 
section we thus used Eqs. (17) and (19). There we will 
demonstrate that the RG data agree quantitatively with 
the DMRG data for U as large as the band width. 



stronger impurity, i.e. smaller p the asymptotic regime 
should also be reached faster since already the non- 
interacting current is closer to a sinusoidal shape. This 
is demonstrated in Fig. 5 for U =1 and p = 0.25. The 
U = transmission for p = 0.25 is |T(fc F )| 2 = 0.22, i.e. 
this case corresponds to an already fairly strong (bare) 
impurity. 



V. RESULTS 

In this section we present our results for the persistent 
current including interaction and a hopping impurity ob- 
tained from DMRG and RG. To demonstrate that the 
Hartree-Fock approximation should not be used in the 
present context we also show results obtained using this 
approximation. We here exclusively consider n = 1/2. 
In the DMRG and RG we have calculated E (4>) for 
< 4> < 7r at 20 different <fri and used the symmetry 
of the data to extend them to the periodicity interval 
— 7r < <f> < ir. From these data we have determined /(</>) 
using numerical differentiation. To make sure that the 
error obtained from this is small we have not only used 
simple centered differences but also applied an approxi- 
mation scheme using Chebychev polynomials. 45 This ex- 
plains why the 4n are not equidistant (see Figs. 3 to 5). 
Up to very small differences both methods give the same 
current. 

In Fig. 3 NI(<f>) is presented for U — 1 and p = 
0.5, which for U = corresponds to a transmission 
\T(k F )| 2 = 0.64, and different N. Additionally the lead- 
ing 1/N behavior at U = [Eq. (7)] is shown as a dashed 
line. As is obvious the interplay of interaction and impu- 
rity leads to a current which vanishes faster than 1/N. 
The DMRG and RG data agree quantitatively, which 
shows that the functional RG provides a very good ap- 
proximation (for a fairly large U) even if the flow of 
the vertex is neglected. The Hartree-Fock approximation 
gives qualitatively wrong results. This shows that corre- 
lation effects not taken into account by Hartree-Fock are 
very important. We thus do not consider this approxi- 
mation any further. For U = 1.733333, i.e. close to the 
phase transition into a charge density wave groundstate, 
the DMRG and RG data are shown in Fig. 4 for the same 
p. For this larger U the suppression is even stronger but 
the RG data still show an excellent agreement with the 
high precision DMRG results. 

As discussed in the introduction in the large N limit 
bosonization predicts that I((f>) is suppressed by a factor 
jq-ciB-t instead of A^ 1 and proportional to sin0. With 
increasing U, as gets larger and one thus expects this 
asymptotic regime to be reached faster. Even though we 
observe a suppression which is stronger than A^ 1 from 
Figs. 3 and 4 it is not clear whether a power-law behavior 
is realized. Furthermore the data even for U = 1.73333 
are still far from being proportional to sin</>. For a 




FIG. 3. Persistent current NI(<j>) for U = 1, p = 0.5, and 
different N. The filled symbols are DMRG data, the open 
symbols RG data, and the solid lines are obtained from a 
Hartree-Fock calculation (for the same TV). The dashed line 
is the 17 = result Eq. (7). 
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FIG. 4. The same as in Fig. 3, but for U = 1.73333, p = 0.5 
and without Hartree-Fock results. 
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FIG. 5. The same as in Fig. 3, but for U = 1, p = 0.25 and 
without Hartree-Fock results. 



behavior of the k — 2 and 3 Fourier coefficients seems to 
be consistent with the above predictions obtained using 
the one-particle language. For p = 0.5 - an impurity of 
intermediate strength - even for large U — 1.73333 and 
up to N — 128 respectively N — 256 lattice sites the 
DMRG and RG data still show a strong curvature. For 
this case the bosonization prediction cannot be demon- 
strated unambiguously although the data show a ten- 
dency towards a behavior which is consistent with this 
prediction. Results for U = 1, k = 1, and different p 
are summarized in Fig. 8. The bare transmission for 
p = 0.9 is \T(k F )| 2 = 0.99 at U = and for p = 0.1 we 
have |T(fci?)| 2 = 0.04. These findings are consistent with 
the results obtained for a local observable - the spectral 
weight close to the impurity. 33 ' 34 



To analyze our data in more detail we have numerically 
determined the coefficients Ik (for k = 1, 2, 3) of a Fourier 
expansion 



I(<f>) =J2 h sin (k<f>) 



(20) 



fe=i 



for different U, p, and N. From bosonization one ex- 
pects, that Nik for k = 1 decays asymptotically as 
N~ aB and that the higher order Fourier coefficients die 
off even faster. To better understand the behavior of 
Nik for k = 2 and 3 we will use the one-particle picture 
mentioned in the introduction. According to this many 
bosonization results can be understood in terms of an 
effective \T(kp)\ which is suppressed by a factor N~ aB 
compared to the non-interacting transmission amplitude. 
Using this and Eq. (8) one concludes that Nik (at least 
for k = 2 and 3) should decay as N~ kaB . This argument 
is based on the use of a single particle picture and needs 
further justification. Fig. 6 shows In [N Ik] {k = 1,2,3) 
as a function of In [TV] for U = 1, p = 0.25 and Fig. 7 
contains data for U = 1.73333, p = 0.5. In addition to 
the DMRG and RG data straight lines with slope —kas 
[see Eq. (3)] (solid lines) and — ka^ G (dashed lines) for 
k = 1,2,3 are shown. Since the flow of the two-particle 
vertex has been neglected in the functional RG, within 
this approach we can only expect exponents to be cor- 
rect to leading order in U. a^ G is the boundary expo- 
nent of the local spectral weight (close to the boundary) 
as we have determined it in Refs. 33 and 34 using the 
RG for much larger system sizes. As expected it agrees 
with oib only to leading order in U. For U = 1 we have 
a| G = 0.291 instead of a B = 1/3 and for U = 1.73333, 
a| G = 0.475 instead of a B = 0.667. From both fig- 
ures it is clear that higher order Fourier coefficient de- 
cay faster than lower ones. On the other hand only the 
data obtained for a already fairly strong bare impurity 
p = 0.25 follow at large N lines with the expected ex- 
ponents (see Fig. 6). Even in this case the DMRG data 
for larger N have to bend further down to reach the pre- 
dicted behavior (slope ~ —kas)- For large N also the 



£ -4 




ln[N] 



FIG. 6. In [NI k ] as a function of In [N] for U = 1, p = 0.25, 
and k = 1,2,3. The filled symbols are DMRG data and the 
open symbols RG data. The solid lines have slope —kas and 
the dashed ones slope —ka% G . The dashed-dotted lines are 
guides to the eyes. Open symbols which are not seen are 
hidden by the filled ones. 




ln[N] 

FIG. 7. The same as in Fig. 6, but for U = 1.73333, p = 0.5. 
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ln[N] 

FIG. 8. The same as in Fig. 6, but for U = 1, k = 1, and 
different p. 

One can stress the single particle picture even further 
and use the non-interacting expression for I(<f) [Eq. (7)] 
with an effective transmission amplitude |T e ff(A;jr)| as a 
parameter to fit our data for different U, p, and N. For 
small ./V this fitting procedure cannot be expected to work 
as good as for larger ones since Eq. (7) only provides the 
large N behavior of the current in the non-interacting 
case. One then expects that |T c ff(&F)| for N — ► oo scales 
as N~ aB (respectively N~ a B for the RG data). Similar 
to Fig. 8 the effective transmission amplitudes show this 
behavior only for small p and large N . For intermediate 
to weak impurities the scaling limit is not reached. We 
have demonstrated that the use of the single particle lan- 
guage can be useful in the present context. Nevertheless 
it involves arguments which are "hand waving" and we 
thus do not follow this route any further. 



as the persistent current. This holds for interactions as 
large as the band width. Our results show that even in 
the cases where the universal scaling limit is not reached 
correlation effects included in the RG are still very im- 
portant. This is demonstrated by a comparison to results 
obtained with the Hartree-Fock approximation which are 
qualitatively wrong. 

Besides the hopping impurities we have also studied 
the case of a single site impurity 33 for selected parameter 
sets and obtained results which are comparable to the 
ones presented above. 

The functional RG can be applied to other microscopic 
models (including e.g. long range interaction), at finite 
temperatures, and for models with disorder. We ex- 
pect that including the flow of the two-particle vertex 

- as discussed for the local spectral weight in Ref. 46 

- will further increase the accuracy of the RG data for 
the persistent current. The flow of the interaction ver- 
tex is of special importance for models with spin since 
in Luttinger liquids with spin the flow of the electron- 
electron backscattering has to be taken into account. 47 
It should thus be included in further applications of the 
RG method on one-dimensional models. Considering cur- 
rents driven by a small difference in voltage within the 
RG it is also possible to include leads and to realize a 
great variety of electronic circuits. It can be expected 
that in small conducting molecular systems of the size of 
a few ten to a few hundred lattice sites which in the near 
future will be accessible to experiments correlation effects 
will be important but the asymptotic limit described by 
the effective field theory will not be reached. The RG 
is an ideal tool to investigate the electronic properties of 
such systems. 



VI. SUMMARY AND PERSPECTIVES 
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